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SUMMARY 


This paper presents an in-depth study of the Thermal 
Response Factor Method enhanced by G. P. Mitalas and J. G. 
Arseneault. A concise analysis of the theoretical mathe- 
matics and physical interpretation of response factors will 
be presented. It will be shown how a new root calculating 
procedure was developed to substantially reduce the computat- 
ional time requirement of response factors by an average of 
88.5%. Comparison testing between the original and the new 
procedures will be presented to validate the new root proc- 
edure. 
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INTRODUCTION 


The Thermal Response Factor Method is recognized by the 
American Society of Heating, Refrigerating and Air-Condi- 
tioning Engineers (ASHRAE) as the current state-of-the-art 
technique to use for the computer simulation of hourly tran- 
sient heat conduction analysis. This method calculates 
hourly heat transfer coefficients (response factors) by using 
the material properties of a particular planar construction. 
These response factors are then combined with an appropriate 
temperature history curve to yield a series solution for the 
heat flux at a given. time. This Response Factor Method is 
the heart of many current energy analysis programs (such as 
NASA's ENERGY COST ANALYSIS PROGRAM - NECAP) . 

This paper is an in-depth study of the Response Factor 
Method enhanced by G. P. Mi talas and J. G. Arseneault, and of 
the RESPONSE FACTOR PROGRAM (RESFAC 1 ) . As this research pro- 
gressed, it became quite apparent that the calculation of 
thermal response factors was being handled in a 'black box' 
manner. To rectify this, a concise analysis of the theoreti- 
cal mathematics and physical interpretation will be presented. 

The second part of this paper will present an investi- 
gation of the response factor's computational time require- 
ment which has been criticized for taking too long. In 
particular, the root finding procedure has been singled out 
as the main culprit. To determine the validity of this 
criticism, the RESFAC code was analyzed and test cases con- 


1 



sisting of various planar constructions were run. No basis 
for the computational time criticism could be found in the 
test cases' results. (See TABLE I) 

However, the analysis of the FORTRAN code revealed 
several portions of inefficient coding. Moreover, it was 
determined that a faster, more efficient, and more accurate 
root calculating procedure could be developed. Upon cleaning 
up the code and installing a new root procedure, the modified 
RES FAC was benchmarked using the test cases mentioned previ- 
ously. The validation testing shows a substantial reduction 
in RESFAC’s conputational time requirement. (see TABLE I) 
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NOMENCLATURE 


c : specific heat (BTU/LB-°P) 

3 

s density of material layer i (LB/FT ) 
f : heat flux function (BTU/HR-FT 2 ) 

F : Laplace transform of heat flux function 
k : conductivity ( BTU/HR-FT- ° f ) 

NRT : number of response factor terms calculated 
P : Laplace transform parameter 

2 

QI : inside surface heat flux (BTU/HR-FT ) 

QO : outside surface heat flux (BTU/HR-FT 2 ) 

R^ s thermal resistance of material layer i (HR-FT -°F/BTU) 

RN : number of real, positive roots calculated 

R : total thermal resistance of all material layers 
( HR-FT 2 -°F/BTU) 

SH^ : specific heat of material layer i (BTU/LB-°F) 
t : time 

TI : inside surface temperature (°F) 

TO : outside surface temperature (°F) 
u : Heaviside unit function 

U : overall heat transfer coefficient ( BTU/HR-FT 2 - °F) 
v : temperature function (°F) 

V : Laplace transform of temperature function 
W : real, positive roots of B(W^)=0.0, k=l, 2, ... NRT 
X, Y, Z : thermal response factor series ( BTU/HR-FT 2 -°F) 

XK^ : conductivity of material layer i ( BTU/IIR-FT- °F) 

XL^ : thickness of material layer i (FT) 
a : diffusivity (FT 2 /HR) 
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n : number c?f material layers 
$ : pulse function 

<j> : Laplace transform of pulse function 
6 : time increment 
p : density (LB/FT 3 ) 

BETA. : reciprocal of the square root of Fourier number 
1 for material layer i 



A, B, C, D: overall conduction matrix elements where 

/A B\ /a, b,\ /a, bA /«„ b n \ 

\c D/ \e, d,y \c» di/ “ ^ Cj1 <•„ / ; 

H = number of material layers; a^, b^, c^, and d^ 
are the conduction matrix elements for each 
material layer i= 1, 2 , . . . , ri 


A 1 



dA, dB, dc, dDi overall derivatives of the conduction 

matrix elements where 




NOTES 




ANALYZING RESPONSE FACTORS 


Physical Interpretation 

Thermal response factors can be defined as time depen- 
dent heat transfer coefficients which depict an object's 
action upon a particular outside surface/inside surface 
temperature difference over a period of time. This action 
is composed of three terms — the X; X, and 2 response factors. 

Using the heat flux eguation 0 “ U *AT along with the 
response factors, the time dependent heat flux values for 
the planar construction shown in Figure 1 can be obtained. 



FIGURE 1 

The outside surface heat flux (QO) at any time t is 

9 

calculated by : 


00 . = 


NRT 

£ 

i=0 


NRT 


( TO . . * X. ) 

t-1 l 7 


-Z 


(TI 


t-i 


V 


i=0 


where NRT is the number of response factor terms calculated, 

TO, . and TI, . are outside surface and inside surface tem- 
t-i t-i 

peratures at times t-i hour. 



QO t = QI fc = U (TO - TI) 

TO = TO t for any time t 

TI = TI fc for any time t 

to be satisfied, the individual summation of the X, Y, and Z 


response factors must equal the overall heat transfer qoef- 
g 

ficient U : 






Mathematical Derivation 

For a homogeneous, isotropic solid whose thermal con- 
ductivity is independent of temperature, the one dimensional 
equation of conduction of heat with no heat generation is: 


( 1 ) 


3 2 v 1 3v 

TT “ *a*Sr “ 


v 

where a* the diffusivity, = — ; and k is the conductivity, 
p is density, and c is specific heat. 

For a planar construction that experiences a steady, 
periodic temperature v(t) , the boundary conditions are: 
v(t) = v 0 <J>(t) for the surface at x = 0 
v(t) = v L <j>(t) for the surface at x = L 
where L is the thickness of the construction. 


(j>(t) is a pulse function of period 26 (6 = discrete 
hourly time interval) and has the values: 


= j , 0 < t < 6 

4> - | ~ | [u(t - 6)] ,6 < t < 26 

♦ * | ~ | C u(t " 6 >3 + j [u(t - 2 6)] ,t > 26 

where u is the Heaviside unit function? 

The heat flux can then be expressed as: 

(3) = 

where R, the resistivity, = i. 

The Laplace transforms of equations (1) , (2) , and (3) 

2 

are now needed. Using the Laplace transform equation 
S = / o e ^S(t)dt, 


n 



where P is the Laplace transform parameter , the solutions 
are: 

,4, s!|-|v.0 


dx 


? = 


(5) 


1 

3p 


,0 < t < 6 


(_ 


w- 


2e 


-P 6 




6 P 


2e- p{ ^ e- 2PS 

T 


6P 


6P 


,6 < t < 26 

,t > 26 


(6) F — i yg 

The solutions (only the real parts) to equations (4) 
and (6) are: 


(7) 


V x <p = V 0 « COSgx f 


F q R SINgx 


_ V <J>g SINgx _ 

F x r + F 0 C0S< 3 X 


where q = / — . This solution can then be expressed in matrix 
•y a 

form: 


( 8 ) 


/v\ / c< 

V F x / V-] 


| SINgx 


! COSgx 

-3 SINgx COSgx 
In terms of heat flux, the matrix can be manipulated into: 



(9) 



SINgx R SINgx 


g q COSgx 

R SINgx R SINg: 



When the planar construction is composed of varying material 
layers, the matrix (8) would be represented as: 
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where 



for i = 1, 2, ... n (number of layers) 


is each layers 2*2 matrix an$ f, is the thickness of the con- 
struction. Note that this solution assume^ that tfysre is 
no contact resistance between the layers. Matrix (9) can 
be represented by: 


( 11 ) 


/V \ / — i\/W \ 

/ Q 1 / B B If VV P I 

\f l / l zl 2 l\W L y 


The solution to matrix (11) is obtained by applying 


the Inversion Theorem for Laplace traps^i^s: 

- ■ A- 52 {£’“ 


In order to simplify this calculation, express the trans- 
formed flux F as a generalized equation; 

< 12 > F - 1 5W 

where <J> is the Laplace transform of the pulss functiop <J>, 


B (P) = “t 


SINgx 


, and M(P) is either 1, A, or D. Bince R t is 


the total resistivity pf tfte piapar construction and U = — 


then B (P) = S uq^* ' 
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To find the inverse of equation (12) , the Calculus 
of Residues* will be utilised: 

1) determine:’ which components in (12) will make the 
equation undefined — 

a. Since If is composed of element (refer equation 
(5)), if P = 0 then (12) becomes undefined 

b. Since (12) is composed of element if 

B (P) = 0 then (12) becomes undefined 

2) dealing with the $ part first — 

a. cancel out the element i in ^ and obtain : 


(13) 


• <$> • = 

•?' = 
•J' = 


1 

6 

1 

X 

1 

X 


,0 < t < 6 


2e 


-P6 


2e -P6 e -2P6 

T + 7 


,6 < t < 26 

f t > 26 


b. By the Calculus of Residues, part I of the 
answer is therefore: 


(14) 


LIM d 
P-*-0 dP 


/e Pt M(P) 

V B (P) / = 


i_ u dM U M(P) dB| 

L UM(p ) + 6 dP - “T— apjp =0 


[ U 2 M(P) dB U dM 1 

. 6 dP " 6 dP J p=0 


,0 < t < 6 


,6 < t < 26 


0, t > 26 
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3) for the B(P) = 0 part — 

3 

a. By using the Heaviside Expansion Theorem, part II 
of the answer is: 



t> 26, j = 3, 4, ..., n= number of hours 
In the above equations, RN is the number of roots of 
B (P) = 0 and W r is the negative of each root (W r = -P). 

When equations (14) and (15) are combined, when 6 is 
set to 1 hour, and when M is replaced by D, 1, and A, and 
evaluated, the equations for the X, Y, and Z (respectfully) 
response factors are obtained; 
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The equations for the X response factors are: 

«U-r • J|fe) • \ • (i)„; '-“"I 

' E“ * ®U.o] * [” 2 • *E[(^) • \ • (f) B • (l- 2 e* eMkl ) . 

-£[(^) * x * (i)„ K - (i - ( ^ ,)3 - 

i=2, 3 , 4...NRT 

The equations for the Y response factors are: 


• -ay 

£f(^r) *(§f) * (i-* ,4Wk) ) 2 * e-C (1+1 > 4w i0j 

k » l L k w j , J 


h “ 


i=2 , 3, 4...NRT 

And the equations for the Z response factors are: 


■ * [» * e»„.J -P • h * (i)„ K • •" tWk ’] 

- f ° * s ».,.«] ' [° 2 • » *£[(^) * % *(!)* (i - ! ” (, “ k,) • •''“"’J 


i=2 , 3, 4...NRT 
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Conduction Matrix Element B 


To begin with, define a purely arbitrary wall composed 
of eight inch thick concrete block (outside surface) , air 
space, gypsum board, and inside air. Using the material 
characteristics described in APPENDIX A, the equation for 
the overall conduction matrix element B will be obtained. 

The equation for B is determined by the following 
matrix algebra: 



( a i ^ / a 2 b 2 \ ^ / a* bj\ ^ /a* b v 

Ci dj / \ c 2 d 2 J \c 3 d s / \ c* 


^ B = (aj*a 2 + b,* c 2 ) * (a s *bi, + b 3 *d„) + 
(a,*b 2 + bj*d 2 ) * ( Cs*b|, + a 3 *d„) 

As generalized equations, each layer's matrix values 
are determined by: 


NON-AIR LAYER 


for W>0.0: 

a i = C0S(/IT* BETA.) 

b. = R i *sin(/JT * BETA.j)/(/[r * BETA^ 

c. = -(/IT* BETA. )*sin(/ir * BETA.)/R. 
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da. = BETA. * sinfv'W'* BETA.) 
2.0 *AT 


db i = ( R J ) * ( 

\2.0*W / \ 


sin(/5T * BETA.. ) cos(/fT* BETA i 


AT * BETA • 


de. ^BETA.^ ^ 


n(/tr* BETA 
/IT 


.) + BETA.*cos(/VT * BETA.)^ 


dd.. = da^ 


as W-+-0.0 

A = 1.0 
B = R t 
C = 0.0 
D = 1.0 

da. = BETA.*BETA./2.0 

l l i 

db i = *BETA . *BETA-/6 . 0 
dc i = BETA i *BETA 1 /R i 
dd^ = da.. 
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AIR LAYER 


a . = 
x 

b. - 


1.0 


R. 

X 


c. 

X 

d. 

X 


= 0.0 


a. 

x 


Upon substituting the appropriate expressions for the 
a's, b's, c's, and d's (keeping in mind that the second and 
fourth layers are air) , the B equation then becomes: 


B = 


^COS ( /W*BETA X ) * [^R 4 *COS (/W*BETA 3 ) + 

r 3 *(/w*beta 3 ) _1 *sin(/w*beta 3 )J J + 

£ £r 2 *cos ( /w*beta 3 ) + 

R 1 *(» / W*BETA 1 )“ 1 *SIN(» / W*BETA :l ) J * 
[-R 4 * (R 3 ) “ 1 */w*beta 3 *sin (/W*BETA 3 ) + 
COS (/W*BETA 3 )J J 
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The only unknown value in the B equation is the variable W. When 
0.0<W<®, B graphs out as: 


B(W) 



FIGURE 2 


In order to calculate the thermal response factors (X, Y, 
and Z) for this planar construction, the values of W which 
makes B(W)=0.0 must be determined. As can be seen in Figure 
2, this complicated wave has no predictable period. There- 
fore, a mathematical root finding technique must be used to 
solve B(W)=0.0. 
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CALCULATING ROOTS 
Original Root Procedure 


The calculation sequence for this method is: 

1. Calcualte first root >30.0 using sequence 3; 
and set ROOT(l) = to this root; set LAST =1 

2. Determine root intervals for locating all roots 
<30.0 but >0.0 

a) if LAST = 1 then set W1 = 0.0001 
else wi = ROOT (LAST - 1) 

b) W2 = ROOT (LAST) 

c) W2 = W2 - 0.0001 

Wl = Wl + 0.0001 

d) DELTA = (W2 - Wl)/20.0 
J = 1 

IC = 0 

e) W3 = Wl + (J* DELTA) 

f) calculate B(W1) and the derivative of B, DB(W1) 

g) calculate B(W3) and the derivative of B, DB(W3) 

h) if B(W1) * B(W3) <0.0 then have a root interval 
so set W2 = W3 and go to 3 

i) if DB (Wl) * DB(W3) <0.0, set IC = IC + 1 

j) if IC - 2 >=0.0 then have a root interval so set 

W2 = W3 and go to 3 

k) J = J + 1 

l) if J >19 then no root interval so set LAST = LAST - 1 
and go to 4 

m) B (Wl) = B(W3) 

DB (Wl) = DB(W3) 

W3 = Wl + ( J*DELTA) 
go to 2g 
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3. Determine the root that is in the interval (Wl, W2) 

a) N = 1 
J = 1 

BP 3 = Wl 

b) calculate B(W1) 

B2 = B (Wl) 

c) DELTA = (W2 - W1)/(20.0*N) 

d) W3 = Wl + ( J*DELTA) 

e) if W3 - W2 >0.0 go to 3i 

f) N = N + 1 

g) if N >25 then there is no root <30.0 so go to 4 

h) set J = 1 and go to 3c 

i) calculate B(W3) 

B3 = B(W3) 

j) - if B (Wl) * B(W3) <0.0 go to 31 

k) J = J + 1 
B2 = B3 
BP 3 = W3 
go to 3d 

l) A = BP3 
B = W3 
FA = B2 
FB = B3 

m) perform Bisection Method on interval (A, B) until root 
is located (i.e. ABS (B (A) - B (B) ) <1.0E-8), 

go to 4 
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4. LAST = LAST + 1 

ROOT (LAST) = ROOT (LAST - 1) 
ROOT (LAST - 1) = W3 

5. Go to 2 

6. End 


Although the Bisection Method guarantees a root once the 
root interval is found, it requires a significant amount of exe- 
cution time to actually locate the root. Compounded with the 
time needed to locate each root interval, this method is quite 
inefficient. 
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New Root Procedure 


The basis for reducing execution time and cost lies in two 
categories: originating a faster method to locate the root in- 

tervals and using a more efficient algorithm for determining 
the roots. 

Root Intervals : 

As was shown , the temperature flux matrix element B 
is mainly composed of sines and cosines. Taken individually, each 
non-air layer's sine and cosine equations have calculatable roots. 
Therefore, this result will be used to segment the B wave (Figure 2) 
into possible root intervals. 

For each non-air layer, use its BETA value to determine the 
segmenter roots (i.e. /IT) . 

RVC = ( i -0. 5 )* x /BETA (for cosine) 

RVS = i* x/BETA (for sine) 

for i = 1 ,2 ,3,4 

Upon doing this for each non-air layer and then sorting all of 
the RVC and RVS values in ascending order, an array V is obtained 
that can be used for root interval checking. 
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^ ~ (Vj, V2* Vj, ••••) ,v,<v 2 <v . 3 

where {V 3 , V 2 } is root interval 1 
<V 2 , V3} is root interval 2 
etc. 

Determining the Roots of B : 

Having determined the root interval checking array, the 
actual roots can then be found. 

The calculation seguence can be summerized as : 


a) set i = 1 

b) set Gl=V i 

c) set G2=V i+1 

d) calculate B(G1) 

e) calculate B(G2) 

f) if ABS(B(G1)) or ABS(B(G2)) ' _< 1.0E-8, a root has been 
found so go to h. 

g) else check the possible root interval at one- tenth segments: 

1) set D = G2-G1, FD = 0.0 

2) FK = FK+1.0, if FK >10.0 then the entire interval 
has been checked so increment i and go to b 

3) else calculate segment end point: 

GNEW = PINT i + ( D*FK/10. 0) 

4) calculate FNEW = B(GNEW) 

5) check sign of FNEW against sign of B { G 1 ) : 

if FNEW*B(G1) >0.-0 then no root in this segment, so 
reset segment front point, G1 = GNEW and B(G1) =. FNEW 
then go to g2 
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6) else there is a root in this segment; store segment, 
front point values: GX = G1 and FX = B(G1), go to h 

h) attempt to locate the root by first using the Secant method; 
if the root 'guess* ever exceeds the end points of the seg- 
ment being checked, then this method failes so use the more 
time consuming Bisection method on the segment end points 
(i.e. GX and GNEW)to locate the root; go to i 

i) the root value obtained was actually the square root of the 
root value so compute the actual root (i.e. ROOT = root value 
*root value); if ROOT >30.0 go to k 

j) reset segment front point: G1 = GNEW and B(.G1) = FNEW; go to 

g2 to check next segment in this interval 

k) end 
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VALIDATION OF NEW ROOT PROCEDURE 


Ten test cases were run to compare the original RESFAC 
against the new program. Input requirements and output 
for both programs are in APPENDIX A. Planar construction 
ranged from a simple, two material layer wall to a complex, 
ten material layer roof. A wide sampling was used for com- 
parison testing and execution timing. The data input and 
both program's output for the test cases are given in APPEN- 
DIX B. 

All testing was done on a high speed. Control Data Corp. 
CYBER 173 computer. Both programs use single precision 
FORTRAN code and were run on a time sharing, batch mode. 

TABLE I lists the execution times and costs for the test 
cases. 
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TABLE I 


EXECUTION TIME & COST COMPARISON* 


TEST 

CASE 

EXECUTION** 

TIME 

% 

REDUCTION 

EXECUTION*** 

COST 

% 

REDUCTION 


OLD 

NEW 



OLD 


NEW 


1 

0.965 

0.228 


76.4 

0.57 


0.45 

21.0 

2 

5.642 

0.131 


97.7 

0.61 


0.45 

26.2 

3 

1.712 

0.292 


82.9 

0.58 


0.45 

22.4 

4 

3.566 

0.415 


88.4 

0.60 


0.45 

25.0 

5 

2.181 

0.344 


84.2 

0.58 


0.45 

22.4 

6 

17.584 

1.341 


92.4 

0.74 


0.46 

37.8 

7 

6.767 

0.696 


89.7 

0.63 


0.45 

28.6 

8 

5.136 

0.533 


89.6 

0.61 


0.45 

26.2 

9 

5.618 

0.601 


89.3 

0.62 


0.45 

27.4 

10 

17.314 

1.034 


94.0 

0.73 


0.46 

40.0 


Average 

reduction 

in 

execution 

time 

was 

88.5% 



Average 

reduction 

in 

execution 

cost 

was 

27. 7% 



* All runs were made on a Control Data Corp. CYBER-173 

** Execution time is in Central Processor (CP) seconds 

*** Execution cost is in System Resource Units (SRU) 
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CONCLUDING REMARKS 

From the analysis performed, it was found that thermal 
response factors are time dependent heat transfer coeffi-r 
cients whose derivations require complex mathematics. In 
deriving the response factor equations, it was necessary to 
calculate the Laplace transforms and the inverse Laplace 
transforms for the heat conduction equation and the heat 
flux equation - both of which are acted upon by a pulse func- 
tion comprised of the Heaviside unit function. Furthermore, 
the keys to obtaining the inverses were found to be the 
Calculus of Residues Theorem and the Heaviside Expansion 
Theorem. 

Testing of the original Response Factor Program re- 
vealed that its computational time requirement was not 
exorbitant. Thus, no basis was found for the criticism 
expounded against it. However, an analysis of the program 
code did reveal inefficiencies in the root calculating 
procedure for conduction matrix element B. This warranted 
the development of a new root calculating procedure. 

From an analysis of the response factor derivation, it 
was decided that each material layers SIN and COSINE com- 
ponents could be utilized in an effective manner. Since 
these components are periodic functions, their zero roots 
are easily calculated. Once these component roots are ob- 
tained for each layer and then ascendingly sorted, a root 
interval checking array can be obtained. Segments from this 
array (values at array positions i and i+1) can then be 
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sectionalized into ten parts and each of these mini-segments 
can be checked to determine if it contains a root for matrix 
element B = 0. This approach to a new root calculating 
procedure proved to be very appropriate. 

In comparison testing, the new Response Factor Program 
executed more efficiently and more accurately than the orig- 
inal RESFAC. The new program's computational time require- 
ment was much less that that of the original program - a 
substantial 88.5% average reduction. 
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NOTES 


25.3 



APPENDIXES 
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APPENDIX A 


PROGRAM INPUT/OUTPUT 

Input to the Response Factor Program (FIN 4.8 extended) 
consists of: 

a) number of material layers 

b) for each material layer (arranged- from outside to 
inside) , its: 

i) thickness (FT)? if air layer then input 0.0 

ii) thermal conductivity (BTU/HR. FT. F) , if air 
layer then input 0.0 

iii) density (LB/FT 3 ), if air layer then input 0.0 

iv) specific heat (BTU/LB.F) ? if air layer then 
input 0.0 

v) thermal resistance (HR. FT 2 . F/BTU) ; if NON 
air layer then input 0.0 

vi) alpha-numeric description (30 characters) 

Output from the Response Factor Program consists of: 

a) echo of the input 

b) thermal conductance of planar construction 

c) response factors series X, Y, and Z 

d) number of hours to reach common ration 

e) number of response factors per set 

f) common ratio 

The common ratio is defined as the common rate at which the 
response factors tend to decrease: 

X i+1 Y i+1 Z i+1 
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APPENDIX B 
TEST CASES 
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t4 O' 


TABLE 1 (Cont.) 


ROOTS 


- hoot # - — old — 

— -0.32755449581283- 

- lr.93392856508856- 

-5.10625948098138- - 

— 9,86100164604460 - 
-16.19980403179670 
-24.123014431-18963 
-33.63074066210538 


0 (ROOT) 


FOR- ROOTS--# OLD 

1 ■ 1 • 191627E-8 

2 7.966970E-9 

3 -2.796240E-9 — 

4 -1.886340E-9 

5 -3.086200E-9 

5.348400E-10 

1.191860E-9 - 
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NEW 


0.32755450010252 

- 1.93392855643586— - 

-5.10625947586985- - 
9.861-00166757376 
16.19980403077625 

24.12301442902663 

33.63074065640694 


NEW - - - - 

-8.96300E-11 

8.20000E-13 

-4.00000E-14 

6.49869E-9 - 

5.00000E-14 

-7.00000E-14 

1.00000E-15 - 



TABLE 2 


DESCRIPTION 'OF CONSTRUCTION' 


,_LAYER ^THICKNESS " CONDUCTIVITY 

_NUfl B ER FT Jl'BTU J»E R I HR M FT I ( F ) 

I .0030 ~ 26.000 

l OlpOOO ~ 0.000'* 


DENSITY ’ ‘SPECIFIC HEAT RESISTANCE 

LB PER CU FT BTU PER ( L B > C F_l (HRUSQ FTMFI PER* BTU 

*80.0 .1000 0.0000 21 SHEET METAL 

0.0 0.0000 .6833 


QT»D — RESF A fl 


* EVpNSE . F ACTOR S 


Mnu» 

i 

Y 

2 

0 

1.60B7382517 

1 .496771 916* 

1. *500055812 

l 

-.2369326768 

•000033638* 

-.0000000063 

2 

0.0000000000 

0.0000000000 

-.0000000000 


Mli-Apo OF MOIJPS REOillPFO TO RflCH* C DM HON*" RATIO • Z 

nu-'r'Er pc OPS PONCE FACTORS ~PEP~~SET '■ ' 3 ~ 

COCHTN' RATIO • .‘0000000000 


NEW 


i55_PP. NSF FACTOR i 


H 0! »R 

X 

Y 

2 .... 

0 

1.66873*251 7 

1 .*58771916* 

1 .*50805581? 

1 

-.239932674P 

.00003365** 

— .0000000063 

2 

o.bboobooooo 

0 . oocooooooo 

-.ooooooooco 


W^>FBJPF__HOU*S FEOUi’PEh TP REACH CONMPN R ATI C_ » 2 

NU * « E R_n f * Vi' P 0 N C E FAC T 0 * $ p E#Jj FT • 3 " 

CQC^PN PAflD ■ ‘ 0. *0000000000' ~ ~ ' ' " ■ 


ROOTS 


OLD — 


-"NO ROOT-IN FALSE" 


"NO 


NEW 

ROOT 00.0" 
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TABLE 3 


0 ESC/ffP TJHH 0 F C 0 NS T RUCT I 


imip 

THICKNESS 

CONDUCT IVITY 

OENSITT 
tft PER CU iFT’ 

SPECIFIC HEAT 

RESISTANCE 



n 

“pTtrpER 'lHP J (FTMFV 

BTU PER (C8HEr“ 

(WCSO FTMFI PER 

B’TO 

i 

.'OP 33 

.*00 

116.0 

• 2000 

0.0000 

STUCCO 

? 

.6670 

.600 

61.0 

.2000 

0.0000 

CONCRETE BLOCK 

3 

“ OT. OOP 6 

' 0.000 

0.0 

0.0000” ' “ 

~ .6050 

INSIDE AIR ~ 


OLD — RESFAC 


P~E SPONS e" FAC TORS 


* oQp t V ' f 


0 

3.313519891*2 

.00*73*90*7 

17026092*003 

1 “ 

-2.0357518186 

.077*69138* 

-.232 5*72853 ” 

2 

— .3181873898 

.122081*7*3 

-.1039*96668 ^ 

3 

-.16732*3708 

.0956263776 

-.06**321910 

A 

-.1057356979 ' 

065 809 0798 

-.0*23797*50 

5 

-.06976*6918 

.0**1981069 

-.0261976595 

6 

-.0*6*561062 

".02 953903 5 0 

'-.0l«e0O32*l 

7 

-.0309938*23 ” 

’.0197221399 

-.0125533739 

A 

-.0206860*98 

.0131650601 

— .008 3790* l 3 

9 

-.013807*819 ’ 

. 00 8 7H7 66 37 

— .005 592909C — 

10 

-.0092163*08 

”.0058657058 

-.003 7332161 " 

11 

-.0061518258 

,00391531 10 

-.002 *918902 


*"»"0FR ps Hnu*5“'«iFOlllREn TO PEACH COMMON RATIO" • 11 

oc prrp'fiNC r~FTCYT5fT">E o *$ft”*'i2 
t ff if r<\— tryt 1 D~ r& St * C 172 C V“' 


NEW RESFA C 


~P~FSPON<c FACTOR'S 


HOUR 


X 

y 

r 

C 


3.3351089055 

.00*73*911? 

1.026892*801 

1 


—2 .0357510095 

.077*691303 

-.2325*72059 

2 


-.3181873091 

.122001*7*3 

-.1039*0666? 

3 


-.16732*3787 

.0956263777 

-.06**821907 

* 


-.105*3*6900 

.0658090799 

-.0*23797*57 

5 


-•0697**6°20 

•0**1981P70 

' -.028197659* 

6 


-.0*6*561*6* 

.029539035 2 

-.01880907*1 " 

7 


-.0309O‘*«*7 5 

“ '.0197221*00 

-.0125533730 ‘ 

-- p 


-.0706*60*99 

.0131650602 

” -.0093790*53 ~ 

9 


— .0130 07*® 

• OOG 70 766 38 

-.0055929090 

10 


‘ -. 00921*3*09 

.0058657058 

-•C03733 21 61 

11 ' 


-.0061518259 

* .0039153119 

-.002*918903 



_NL?*?Er_OF v OtRS FFOIIIPFO'TO PEACH CO^Hn^ ratio"* ii 
NIPPER OF" CFSPONCF 'FACTOkS Vffp s'E T'V 12"" 7 ' 

“CCPPCK RATTP '.667A0172IA 
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TABLE 



FOR ROOTS # OLD 

1 -4.60508E-9 

2 3 • 28262E-9 — 

3 4.50591E-9 

- — - 4 5.74390E-10 

5 3.98570E-10 - 

8 2.42400E-10 

7 5.67000E-1 1 
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( Con t . ) 



TABLE 4 


DESCRIPTION OF CONI TRU CTTON 


LAYER 

-rmmsr 

CofTouCflvlf? 

transit f 

— rp'RrnrTirrr — 

RESISTANCE 


NUMBER' 

FT 

BTU PEP (HR )( FT) CF ) 

LB PER CU FT 

BTU PER (LB) (F> 

“"(HR)! so FT ) ( F ) PER 

BTU ’ 
"CONCRETE" 

1 

1.0000 

.*70 

61.0 

• 2000 

0.0000 

2 

O.CuCO 

0.000 

0.0 

0.0000 

.9100 

AIR SPACE 

3 

.0623 

“*.*20 " 

100.0”“ 

• 2000 

0.0000 

’ GYPSUM BOARD' 

4. 

6'. CO CO 

0.000 

0.0 

0.0000 

.66 30 

" INSIDE AIR ‘ 


OLD RESFAC 


RESPONS E f A CT O RS 


HOUR 

X 

r 

Z 

0 

2.7019910060 

.00000*3*9* 

.902009*93? 

1 

-1.5027901353 

.00137171*1 

-.3592050606 

2 

-.260*62322* 

.010**98*96 

-.09591*1188 

3 

-.1353*16110 

.021796858* 

-.0*256*6*23 

* 

-.087892*779 

.026929*928 

-.02070*6*51 

5 

-.06*322*505 

.0269628958 

-.01969709*7 

6 

-.050**62651 

.02*6772257 

-.015673*50* 

7 

-.0*11015290 

.021690039* 

-.0130731950 

8 

-.03*1619995 

.0107069269 

-.010929*008 

9 

-.020661*781 

.01599*0218 

-.0092032*65 

10 

-.02*200679* 

•01361*6009 

-.0077768730 

11 

-.020*699917 

.0115652090 

-.0063027591 

12 

-.0173351*51 

.00961**69* 

-.0055766122’ 

13 

■“-.01*6089*11 

.00632*7083 

00*7261331 

1* 

-.012*502106 

.007059*1*7 

— .0090061772 " 

IS 

-.01055*1392 

.0059667*98 

- .00339620*7 

16 

"-.0089*7*2*9 ~ 

.0050750956" 

-.0020792386 

17 

* - .0075b 55 566 

’* .00*3028691 ~ 

-.002**10192 

18 

-.006*310775 r 

“‘.0036*80965 

-.002069519* 

19 

-.005*523*55 " 

.0030929*17 

-.00175*5676 





NUMBER Oni'QURS REQUIRED "TO’ RE AC H COMMON RATIO • 19“ 
NUMBER OF RESPGNCE FACTORS PER SET ■ 20 

— - 



COMMGN RATIO - .0*7618302* 
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TABLE 4 (Cont. ) 


NEW RESFAC 


»yspoN T f"y«C To»s 


Htru® t *■ 'v ■ z‘ 


“HT 

2.70l09l03'4l 

TD000045273 ‘ 

.9028094932 


1 

=!• 3827001625 

.0013717363 

-.3392830790 


* 

'-.2604623713 

*.0104498499 

-.0959141192 


3 

-.ruumrt 

.0717468587 

-.0475646477 


4 

-.057*0247«3 

.0269794930 

-.0267646454 


5 

”-.G4432?4«8* 

. 026 962 8958 

-.019807*949 


6 

-.0504462653 

.0246772257 

— • C 159 734506 


7 

-.0411013706 

.0216000374 

-.013073195? 


p 

-.034161 909 5 ’ 

• 0 1 8 7CC 9?69 

-.0109294009 


, Q 1 ■■ 

— .02 865147 31 

.0159940219 

-.0092032465 


1C 

-.0242006776 

•0136146008 

-.0077768730 


11 

"-.020460091 6" ~ 

”.0115652089 

- — ' ’ -.0065827591 


IZ 

-.0173351450" 

" .0098144694 

-.0055766122 


13 

' -.0146**9410 

• OOP 374 70P2 ' 

— .004 7261 53 1 


14 

-. 0124*07104 

. 00705 94146 

-.004006177? 


15 

-;Cl033411Ol 

.005965 7497“ 

-• 0C3 3°62047 


16 

-;0369474?4»» " 

' " .0050750955 

-.0028797385 


' It 

’-;0075?55567 

.0043028690 

“'-.0074410192 


l*' 

-.0064310774 

' .0036480965 

-.0020693194 


19 

’-'.0054573455 

.0030929417 

-.0017543676 








SFOUrRcCTTtrfiglCH coupon -^Tirr. 
“ NU85E*~!jr"5V$P0NeE‘ FACTORS’ »E9 SET -20 

- 

r — 

CCUMGK frATlO • .647 PI *30 11 
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TABLE 4 (Cont.) 


ROOTS 


ROOT# OLD 

1 0.-16508893214859 

2 0.88663631632084 

3- L.82H-T532969845 

4 2 .XO 89X699922X66 

5 4.87369X364X50 1 2 

.6 7.89216455071565 - 

7 1-U6850XX89 151-147 

8 16.2425649X751424 

9 21.56196345126659 

10- 27.64206980186293 

11 34.4817X223186940 - - 


8 (ROOTS) 


FOR— ROOTS- # OLD— 


1 2. 128693E*8 — 

2 l.Z63430E=9 

3 -3.8894Z0E=9 . 

4 1.059484E-8 

5 2.079239E-?8 

6 »1.933367E»8 - - 

X -^1.325572Ew8 

8 2.175127E»8- 

9 1.086137E^8- _ 

*- 10 -s*2.7l 1985E®8 

11 1-.036738E-8-- 



NEW — 

0.16508893366130 

0.8866363157-1600 

1.8211 153271-2176- 

2.70897699132630 

4.87369737239763 

7.89216455605310 

11.68507788855550- 

16.24256497291390 

21.56196345365345 

27.64206980859319 

34.48177223524772 


-NEW 


-3.80000E-13 

— — 7 .OOOOOE-14 

..... ^6*4liaOEwlO ... 
-2.67600E-11 ... - 

-- -7. 00000E-14 
-1.03077E-9- - 

-2.71100E-11 
7.33800E-11 - -- 
-1.09000E-12 
-1.61000E-12 
-4.57000E-10 



TABLE 5 


61 scRimoN of const^c™ 


lAVFR 

TMKKNFSS 

CONDUCTIVITY 

DENSITY 

SPECIFIC HEAT 

RESISTANCE 


KH**Pf 0 

FT 

TTu »F> (HR ) ( FTITFI 

LB PE* CU FT 

BTU PFR (L8>m (HPlTSO F f ) ( F ) PER 

BTO 

1 

• 04 1 f 

.830 

55*0 

.4000 

0*0000 

STONE 

7 

• 0?13 

• 110 

70.0 

,4000 

o.ooco 

FELT 

3 

71670'*“ 

• 025 

2.0 

.2000 

O.OOCO* ‘ 

" C E l L U L*A R~*G Ci S $ * 

'4 

•'CO 96 ’ 

26*000 

480.0 

.1000 ~ 

"O.OOCO 

ShEET’HETaC 

9 

0.0606' 

0.000 

0.0 

0.0000 

.6950 

INSIDE AIR 


OLD RESFAC 




HOUR 

t 

V 

2 

6 

I”. 80741 34 526 

•074316 3l 91 

•3484981684 

1 

-1.7660938687 “ 

~ .0549331626 

— .2 17 99 00359 " 

2 

-.0014456604 

.0006165926 

-.0006354967 

3 

-.0000034153 

.0000024075 

-.ooocoSoeol ' 

4 

-.0000000105 

.0000000084 

-.0000000070 

9 

-.0000000000 

,0000000000 

“-.0000000000 


\iiH'aFo n'p HOURS ® Foil f ® f”o to rfach c n h h r ) n ratio ■ 5* 

► * mt* p f r>p p cue F_ fS rroR_s p e £j$ f t « 6~" “ 

>“n>»pn paTIO « .0033832881 ~ * ~ 


NEW RESFAC 


^fVponTe FACTORS' 



~ HOUR 

* 

'Y" 

I 


0 

1.P974 134525 

.0743163181 

.3484891084 


1 

-1.766093*686 

.0545331620 

-.2179600359 


2 

' -.001445^654 ’ "■ 

• OOCfcl 65926 * 

-.0006354967 


3 

-.0000034153 

.0OCO024O75 

-.C000020801 


A . - 

-.000000010* 

".0000000094 

-.0000000670" 


* 

■ -.0000000000 ' ■' '■ " 

.0000000000 

-.0000000000" 



* KU M * a = R 0^ CUI RED TO FEAC'h'CO^HPK' RATIOV '5 

“NTJ^e.rF ~CF _ Fr FP0NCE“F4CTCRS “PEP"SET «~"*6 — 

~'CO*rrN~0ATtO“‘* *00339 3 29 61 
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ROOl 1 -# 


OLD 


5 (Cont.) 


ROOTS 


NEW 



5*68890723527994- 

7.23930076296122 

23*82467391296810 

46*56416592606750 


B(ROOTS) 

- -NEW - 


4.42000E-12 

-3.36164E-9- 

-9.11880E-10 - 

6.80000E-13 
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TABLE 6 


... 



DESCRIPTION of CONSTRUCTION 

.. 

• -• ■ — — 

l A*FP 

TMKKNfSS 

cnsniiPTi vitv 

DENSITY 

SPECIFIC HEAT 

RESISTANCE 


M'«**f* " 

PT ATU (HP M FT) |F) 

LA PE» CU FT 

ATU PER (L A) I F ) 

(H?) CSO FT) tF| PER ATU 

" f " r 

.0400 

.400 

116.0 

.2000 

0.0000 

STUCCO 

7 

• 0630 

.015 

74.0 

.2000 

0.0000 

INSULATION 

3 

.,0400 

.400 

116.0 

• 2000 

0.0000 

STUCCO 

4 

1*0000 

.380 

106.0 

.2000 

0.0000 

NASONARY 

5 

• i 250 

.250 

5.7 

• 2000 

o.oooo 

rigid INSTiLATnJff - 

6 

• 1250 

• 420 

100.0 

.2000 

0.0000 

GYPSUH 

7 

0.0000 

0.000 

0.0 

0.0000 

.6850 

INSIDE AIR 


OLD RES FAC 


RESrONSE_>>CTORS 



HOUR 

X 

”Y " " 

2 


0 

1.4133452868 

-.0000016845 

1.0234028097* 


1 

-1.1675786059 

•0000016986 

-.2922661790 


2 

-.0287241362 

.0000006880 

-.1269399552 


3 

-.0094665212 

.0000158459 

-.0696275927 


4 

-.0058867493 

.0U00965331 

-.0445554251 


5 

-.0045961430 

.04.02893197 

-.0318517154 


6 

-.0038760191 

•0Q05614340 

-.0245613444 


7 

-.oo33i>eee5o 

.0009223559 

-.0199316112 


8 

-.0030234883 

•0012630116 

-.0167516158 


9 

-.0027478187 

.0015717229 

01443*334 4 


10 

-.0025216815 

.0010335561 

-.0126736599 


11 

-.0023349639 ’* 

.0020448307' 

-.0112905141 


12 

-.0021778771 

.0022C80835 

-.0101785102 


13 

-.0020437032 

.0023286077 

-.0092681270 


14 

-.0019276198 

.0024128978 

-.0085119026 


15 

-.0018260430 

.0024660386 

-.0078760512 


16 

-.0017362357 

.0024960470 

**"— • 00733564 j 3 


17 

-.0016560803 

.00250531 51 

-.0068717849 


18 

-.0015936905 

.0024986835 

-.0064698635 


19 

-.0015163328 

.0024795039 

-.0^61164099 


20 

-.0014563340 

.00245052 70 

-.0058083259 


21 

-.0014030280 

.0024139934 

-.0055323309 


22 ~ 

-.0013517122 

.00237171 87 

-.C052646036 * 


23 

-.0013038151 

.0023251670 

-.0050604034 


24 

-.0012508698 

“ ‘ .0022755169 ' 

-.0048559191 


25 

-.0012164939 

.0022237162* “ 

-.0046690514 * 


26 

-.0011763732 

.0021705143 

“ -.0044942787 


27 

-.0011382486 

.0021165137 

-.0043325451 


28 

-.0011019055 

.0020621967 

- .004 1 fa 11 707 


29 

-.0010671655 

.0020079211 ~ 

-.0040367794 


30 

-.0010338796 

.0019539995 " 

-.0039042409 


31 

- • 00 1C01 9227 

• 00 19C 065 53 

"-.0037766241 


'"3 2 

-.0009711891 

»00l64d06 52 

-.0036551595 


33 

-.0009415890 

.0017963638 

* - .0035392uo2 


34 

-.0009130455 

.0017416515 

“ -.0034282382 


35 

-.0008854922 

.0016960017 

-.0033219035 

4 

”3 6“ 

-.0006566716 

.0016474659 

- • C03 2 1952 8 7 " * 

« “ 

37 

-. 00C8331 331 

.0016000783 

- . 0u3 1 2 10956 

i 

36 

-.0008062321 

.00155385 92 

-.0030262333 


39 

-.0007841290 

.0015086160 

- . C02 9 3 4708 2 * 

i 

40 

- • Q0076C78 60 

.0014649553 

-.0026463187 


41 "** 

-.0007381769 

.0014222651 

— *002 7608d89 


39 




TABLE 6 (Cont.) 


sz 

-.00071 62661 

•00136073 3$ 

-.0026782639 

A3 

-.0006950287 

• 00 13 A033 23 

-.0025963063 

A A 

-.000674A394 

*" .0013010936 

- .002 320t>93 A 

A3 

-.00063AA731 

\0012629448 

-.002AA391A2 

A 6 

-.00063311 AO 

*.0012258772 

- .002 3 7 32683 

A 7 

" — .0006163 3 5 3 

•0011698689 

- .002 3026 63 A 

A 8 

-.0003981199 

•00113A6952 

- • 002 2 3 A61 A 5 

A9 

-.OOO50OAA92 

.0011209307 

-.0021694420 

50 

-.0003633037 

• 00 10 d 794 99 

- . 002 10 A2 7A 7 

31 

-.0003A66728 

".0010339273 

- . 002 0 A 20 A1 1 

32 

— .0005303 3AA 

•00102463 73 

*-.0019816770 

53 

-.00051 A9732 

.G0099A63A7 

'-.0019231203 

3 A 

" -.0004996904 ' 

~ .00096535 A3 

-.0016663132 


NUHbt* OF HOURS RE OU IRE D TO~ REACH COMMON RATIO '■ 54* 

HUMBER OF RE SR ON Ct FAC TORS PER SET • 55 

COMMUNJUTIO j« .9703090791 _ 


NEW RESFAC 


— 


*"* rfsponV 

FACTORS 













HOU® 

X 

r 

r 


0 

1 • Al 33A5A20A 

.0000000000 

1 . 023A0161 19 


1* 

— 1 • 1 675 7P961 7 

•0000000016 

-.2922650222 


2 

— .02fl7?Al 355 

. 00000068*4 

1269399500 


3 

-.C09APP51 99 

. 0CC015P‘53 

-. 069627SPPA 


A 

-.005 fl P67A««9 

. 00C09 6 5 3? A 

-.CAA555A213 


« 

- • 004 A 96 1 A 1 P 

,0902593190 

••031 *517119 


6 

-.003" 760 1 A 1 

.000581 A33A 

-.C2A6613A1A 


7 

-.oo3App»»f a: 

.0009273551 

-.C1993160PA 


0 

-.00302PAP7A 

.0012630113 

-•01675161 33 


9 

-.0027A7P1 79 

.001 5717225 

-.01AA353321 


10 

-• 002521 6*07 

. CO 1 83 3 5557 

-.012673656P 


11 

-• 0023 3 A96 3 1 

. 0020A Afl 303 

-.0112905121 


1? 

-.0021 77*763 

" • 00220P0P32 

-.01017*5164 


1 3 

-.0020437025 

.00237P6P7A 

-.0(926* 12 53 


\*> 

-.00192761 71 - 

.002A128975 

— .0^*5119010 


15 

-.001*7P0A2A 

. 00240*0 3P 3 

-.007e7AOA97 

- - .. 

1 6 

— - _ #00Z7 34,71 « 0 

,OO2A96C*60 

-.CC733A0A1P 


17 

-.0016560797 

.00250531*9 

-.006*71 7*35 


10 

-.0015038900 

• 00? A9 06 P33 

—.0^64690622 


* 19 

" *"* — . 00151 a 337? 

• 00? A 7950 3 A 

-.0C611PA0PP 


20 

... - -.OO1A503335 

• 0P?4505?69 

-•C n 50CA37AP 


21 

C01A030275 

• 002A1 3993* 

-.005532337 p 


“*?? 

. -.001 3 51 71 1 P “ 

' .002371 71*6 

-.rC52*A602fi 


23 

-.001303*14* 

.0023251670 

-.C05080C025 


2<. * 

-.00125PP693 

.0022755 1P9 

-.C0AR5591 83 


2 5 

-.00171^4915 

.0022237163 

-.00466*050* 


?l 

-.oon 703775 

.00217051A3 

-.00AA9A2779 


27 

-.00113824*? 

.0021165130 

-.004 332 54 A* 


?P 

-.0011019051 

.0020621900 

-.0041011700 


- ?0 * 

-.0010671*51 “ 

.0020079212 

-.004039778P 


' 30 

-.001033R793 

.0019539996 

-.0C3O0A2A0? 


' 31* ‘ 

— • 001001 0??4 

.CO190C655A 

-.0C3776673P 


3 ? ~ 

-•CCC97tl "AP 

"* .OO10APCE53 

- • OC 2 4 55 1 * 9C 


~ 33 

-.COOQA] 5 0 a 7 

“ .0017963639 

-*003539?07P 


3A 

" "* ‘ ' -.0009130452 

.0017A56517 

-,0C3A?A?37P 


35 

-.OOOPP549I9 "■ 

.0016960019 * 

-.003 3 21 P03 1 


36 

" -.cooe 5PP71 3 

.0016A7A6M 

-.0C321952P4 


37” 

-.000*331 3?A “* 

.001 60007P5 ’* 

-.0C31210955 


”30 

-.000*0P?3!9 

• CO 1 55 30 594 

* -.0030262330 


39"'**' 

-.POC^PM? AA 

" .00150PP1P7 *’* 

"-.0029347079 


AC 

" ““ -.0007A07P75 

.001A6<.9555 

-.CC?P*531 05 

l 

A 1 

— . 0*3073 5 1767 “ * 

*001 4222653 

-.002760**87 
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TABLE 6 (Cont.) 


47 


-.0007162630 


.0013007361 

-.0026782637 

43 


-.0004050295 


•0013403525 

-.0025933062 

44 


-.0004744103 


.0013010053 * 

-.0025208032 

45 


-.0006544750 


.0012620430 

-.0024450141 

4fr 


-.0006351119 


.0012258774 

-.0073732682 

47 


-.000616335? 


•0011*09601 

-.007302*633 

4 4 


-.00010*1199 


.0011549954 

-.0022744144 

49 


'-•0005904491 


.0011209300 

.0021634478 

50 


-.0005633056 


.0010879501 

-.0021042747 

51 


-.0005466727 


".0010550275 

-.0020420411 

5? 


— ^ooosiosm 


"•001024 8 375 

-.0019816770 

53 


—.0005 14*751 


.00C9946549 

-.0019231205 

54— 


- .0004006*03 


.0000653547 

-.0018663132 


MJPBE* OB HOURS REOUIRFO TJTV ■;'* CMC OM HPN _ R R T~j 6 • *'5* 

Tiu^FfofBTSRona 'ritTDRSTrr set • 55 
~CQKfotrffrrr<r= ;.7oso<> 6»'fcT> 


ROOTS 


— ROOT-# OLD NEW 


1 0,029933697 02000 0 . 02993366908454 - 

8 0,21242985580668 - - - 0.21242984890035 

3 0.61496372900272 - - 0.61496373599719 

4 1,09618246581550 - 1 . 09616246155674 - 

5 1 , 66638739781747 - - - — 1.66638739147555 

6 2 , 32695013035205 - •— - 2 . 32695012724621 - 

7 2.96406703559819 - 2.96406703639735 

-8 4 *26934 12077 1537 - 4.26934121302924 

9 5.96470442943337 5.96470443677839 

1-0 7.90782511376565 7.90782511203173 

14 9.38101582073091 - 9.38101581368824 

12 10.74572406709456 — - - 10,74572406078107 

1-3 13.20137627111041 — 13,20137627948367 

14 15.14241773577010 15.14241774014755 

15 16.58063080041711 16.58063079648036 

-16 19.19008477162163 19.19008477625005 

17 20.55331112918213 - 20.55331113165494 

18 23.66926094969699 — - - 23.66926094966527 

19 27.32398281315591 - - 27.32398280793188 

28 29.01735456492725 29 . 01735456421909 - 

21 32.16377783473581 32.16377783283110 



TABLE 6 (Cont.) 


B (ROOTS) 


_..FOR-R00TS n - OLD - NEW 


1 -1.46805208E-6 

2 4.47713380E-T - 

3 1.8824720OE-7 

4 5.84924700E-8 

5 »5.49342000E«8 

6 -2.1l869900E"8 .... 

7 7.39351000E-9 

8 -9.94726500E-8 

9 1 *3502Q560E“7 

— 10 1. 477681 00E-8 

-- -U -3.02108600E-8 

12 -- -2.44352100E-8 — 

13 2.34856100E-8 

14 -5.98533000E-9 - 

15 4.85016000E-9 

16 -4.95923000E-9 

1-7- — 3.20458000E-9 

18 2.62937000E-9 - 


19 -8.93784000E-9- 

20 4.14600000E-9 

21 -6,36227000E-9 - - 


1.00000E-14 

- -9.00000E-14 — 

-1.40000E-13— - 

- 2.24335E-9 

-1.27855E-9 

-1.40000E-13— - 

3.80000E-13 

— -1 •34000E-12 

-3.19200E-10 

4.58000E-12 

- 7.20000E-13 

-4.60000E-13 -- 

-2.11340E-10 

8«OOOOOE*»14 — 
7.00000E-13 
1.48200E-11 - 

1.74200E-10 
2.54531E-92 - - 

- -8.12000E-12 

-8.40000E-12 

-2.40000E-13 
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TABLE 7 


DESCRIPTION ..OFjCdNSTRUCtTON 


LATER THICKNESS CONDUCTIVITY ‘ “ * 'DENSITY SPECIFIC HEAT ’ ASSISTANCE * __ 

NUMBER FT BTU PER (HR ) ( FTMFf “ LB PER CU FT " BTU PER CLBMFl (HRMSO FTXF) PER BTU 

* 1 • 3330 .770 125.0 .2200 ' ‘ ' 0.0000 FACE BRICK 

2 'O.OOCO “ 0.000 0.0 * 0.0000 " .9100' m " ' ” AIR SPACE 

’ 3 * ,0633 " ‘ .070 ' 37.0 .6000 ' O.CGOO * SHEATHING BOARD 

4 .3330 ” .025 ' 5.7 _ *.2C00 O.OCOO INSULATION 

— ^ .0625 ".420 ”'100.0 • 20 CO* '0.0000 ' GYPSUM BOARD’ 

— 6 O.OOCO 0.000 0.0 O.OOCO .6850 " * INSIDE AIR 


6 O.OOCO 0.000 0.0 O.OOCO .6850 " * INSIDE AIR ’ 



OLD RESFAC 



response'fXctors 



HOUR ' X * Y 2 

0 * _ T 5.1759961910 .0000002162' i 8567490 859 " 

" 1 -3.2693245014 " .0001735606 "-.5063325017' 


2 -.7627360815 . 0019 B 7590 e -.1769113994 



3 - 

-.3954482090 

'.0053010003' 

" " ' ' -.0661057275 



4 

- .22 36660060 

.0077112856 

-.0250009737 



5 

-.1372336219 

.0082363141 

-.0106121783 



6 

-.0902026585 

.0075773934 

-.0047067431 



7 

-.0624053962 

.0064230232 

-.0023146240 



6 

-.0446903909 

• OU 51 9409 87 

— .0012 7b536 6 



9 

-.0327070454 

.0040683314 

-.0007067634 



10 

-.0242524535 

.0031605263 

-.0005255624 



11 

-.0181221793 * 

.0024210054 

-.0003 703670 



12 

-.0136018523 ’ “ • 

.0010436046 

* "-.0002693425 



13 

-.0102351670 ‘ 

.0013969911 

~ -.0001993224 



14 

-.0077130398 

.0010593995 

-.0001486973 



15 

-.0058172472 " ‘ " " 

.0000012596* 

-.0001117658 



16 

-.0043895017' 

• 00G6O55652 

-.0000641465 



17 

■-.0033130658 

.0004575046 

' ”-.0000634302 



10 

-.0025009878 " ' " r ' 

.0003455492 

-.0000478500' 



19 

-.0018881264 

.0002609535 

-.0000361113’ 



20 

-.0014255162 

• Ou 0 1 9 705 21 

-.0000272583 



21 

-.0010762808 

• 00014 b 79 16 

“ -;0000205781 



22 

-.0006126172 ~ 

.00011234 76" 

-.C0001 85360 



23 

-.0006135507 “ 

• 00 OOd <t82 87 

- .0000117290 



24 

-.0004632519 " 

• 00006404 90 

- • COOOG00 56 3 



25 

-.0003497721 ’ 

.0000483605 

* ’ ■ ™ ” -.0000066067 



NUMBER OF HOURS RE OUIRED TO REACH COMMON RATI0_«_Z5_ 
NUMBER' OF "RESPCNCE FACTORS PER SET • 26 ' “ ' 
COMMON RATIO • .7550396065 
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TABLE 7 (Cont.) 


NEW RESFAC 


~yfTPPNSg FAC TORS 



HOU® 

V 

Y 

2 


0 

5.1730967394 

.0000002177 

.8567491015 


1 

- 3.2691245749 

.0001735577 

-.5063376014 


2 

-.7627360529 

.0019875912 

-.1769113981 


3 

-.3954482093 

• 005381 8 C 07 

-.0661057268 


4 

-.2236660857 

' .0077112859 

-.0258009732 


5 

-.1377336213 

.0082363143 

-.0106171779 


6 

-. 0902076*78 

. 0 C 75773935 

-. 0 C 47 C 87428 


7 

' -.0624053074 

"“.0064230233 

* -. 0023146237 * 


e 

-.0446903981 

.0051940988 

-• 0 C 12785364 


9 

-.0327070448 

.0040853314 

-.0007867433 


^10 

-.0242574529 

' .0031605263 

■ “ " "—. 0 CO 525562 3 


n 

-.0181221788 

.0024210053 

-.0003703670 


l? 

-. 013 * 018519 * 

.0018436045 

-.0002693425 


..... 13 .. . ... 

-". 0102151*66 

. .0013989911 

-.0001993223 


14 

-.0077130395 

.0010593995 

— • C 001 48 B 073 


15 

—.0059 172469 

.0008012595 

-• 0 C 01 1 1 7 * 57 


16 

— . 0043 P 9501 5 ~ 

~ . 0006055851 ’ 

-.0000841464 


'17 

00331 306 56 " 

* ~ .0004575046 

— . 00006^4302 


18 

-. 002500987 ?* 

.0003455497 

" ” — • CC 004 795 CO * ■" 


19 

-.0018881263 

" . 0002609535 

- -* 0 CC 0361 1 1 3 “““ 


30 

-.0014755167 

.0001970521 

— .00002725 8 3 


21 

-. 0010762 P 07 

• 00 C 14 87916 

-. 0 CC 02057 P 1 

1 

2 ? 

" -. 0008176172 " 

1 .0001123476 

““ *" -. 0000155360 " 

1 

23 

-. 0006135*07 r 

.0000840287 " 

'-.0000117298 

i 

24 

-". 0004*37518 

.0000640498 

— * 0 C 00098563 


- 25 “ 

-.0003497721 

" . 000 C 4 e 3605 " 

— « 0 C 00066 ® 67 


r s~Rrtcrr a rtrrirTTACK c o* pom ratio ~t$ 

NtJ*fBF Rl'O E < f 0 NO r F A £T 0 k % pTp S F T • ‘ 26* 

CC>K0n R AT 3T0 ■ ;735rt9fc& T7 
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TABLE 7 (Cont.) 


ROOTS 


HOOT # 

OLD 

NEW 

l — 

0.20090507216704 

0.20090507502694 - 

2 

0.84220454393237 

0.84228454471546 

3 

1.02137107973186 

1.02137107926540 

4 

2.17346651520421 

2.17346650715290 

5 

5.11260893999551 

5.11260894227845 

6 

- 6.05598433406738 

6.05598433479236 

7 — 

8.22009552003902 . 

- 8.2200^552583131 - 

0 

15.68213748921056 

15.68213749230597 

9 

- - -- 16.65349741351815 

16.65349740592023 

10 

20.10396604661435 

20.10396604274922 

11 

30.63091142382473 

30.63091143841234 


B(ROOTS) 


FOR 

ROOTS # OLD 

NEW 

1 

6 . 7804 32E-8 

2.072OOE-11 

2 

6.7907OOE-1O 

-2.00281E-9 

3 

■7.31 1400E-1 0 

3.30000E-13 

4 

5.621767E-8 

-7.90000E-13 

5 

8.845090E-9 

-4.60900E-11 

6 

=3.960310E-9 

-1.52034E-9 

7 

5.605072E-B 

4.30000E*»13 

8 

t9.308060E~9 

-1 .15000E-12 

9 

-1.537046E-8 

3.25792E-9 

10 

3.210414E-8 

“1 .76000E-12 

11 

5.50O3OOE-1O 

-9.69469E-9 





TABLE 8 


MsJc'ftjyfi on nfjo n s t m t i on ‘ 


AVER THICKNESS CONDUCT I VI TT OEMS TTY SPECIFTc HEAT _ RESISTANCE 

• "••Ft ^ FT' Rtu' PER iHSICFlTTFr _ LB PER CU~FT_JBTU PER CLBMFI* (HR) (SO FTMFI PER BTU" 

_J_ 2 >0100 * ' 2.300 ’ 70, 0' ^ • 3500 " 0.0000 __ ASPHALT SHINGLE 

7 .0A2OJ .067 • 2900 J/O.OOOO PLYWOOD SHEATHING, 

3* .2300 .025 ".5“ *‘.1600 “* 0.0000 INSULATION 

“A 6.0000 0.000 " 0.0 0.0000 ** ■” .9100 AIR SPACE * 

~ 3 , OA 2*0 '.067 3 A .0 .2900** O.OOOO' ‘PLYwOOO 

* 6 .5000 ' “ .025* * .5 .1600 " ' * ~ "0.0000*" IN SUL AT ION * 

7 .0625' . A 20 _ 100.0 ".2000 " 0.0000 *" GYPS'JM BOARD 

p o.oooo o.ooo o.o^ o.oooo'*** .6850' * Ins ioe air 


OLD RES FAC 


Ri5ZQN5.ElLA£JflRSl 


HOUR 


X 

Y 

1 

0 


.7251569790 

.0003756736 

•037 AA 57518 

l 


-.6562617879 

•0036660986 

-.5326726953 

7 


-.0101950075 

.0055020236 

— . 175 913762 5 

3 


— . 00 7A7 98 669 

.0050052262 

-.0596083780 

A 


— • 005507 2 1 A A 

.0039969156 

-.021 1351227 

5 


-.00A05653A3 

.0030A53759 

-.008 2261789 

6 


— • 00299 fl 5A 57 

.0022765117 

-.003 663A230 

7 


-• 00 22019 1 5A 

• 00168 798 8 9 

-.0019206027 

B 


-.0016223972 

• 001 2A 7 20 1 B 

-.001 1622678 

9 


— .00 1 195 A 206 

.0009200953 

-.000 77A268 2 

10 


-.0008308205 

.0006783189 

-• 000 5A383A 6 

11 


• 0006A 90 1 60 

• 000 A99 9? 5 A 

-.0003920535 

12 


-.OOOA76215?" 

.0003683998 

-.0002860650 

13 


00035?365i 

000271A615 

— .000 209B68 6 

1A 


0002 606 3 A 1 

.0002000258 

- .000 15 A 3 A 1 A ’** 

15 


-.0001913070 

• 0001 A 7 3869 

-.0001136275 

16 


— • 0001A09613 

. 0001 0 R6000 

-.0000836932* 

17 


- . 0001 03 B 6A9 

•0000800202 

~ • 000 061657 8 * 

IB 


-.0000765311 

.0000589616 

” 000 0A5A 28 2 * 

19 


-*.0000563906 

.0000A3AAA8 

— • 000 03 3A 71 9 

20 


-.0000615605“ 

.0000320116 

— .000 02 A 662 9 

21 


-.0000306158' 

.0000235072 

-.0000181723 


NUHVf S R E QUI P E D TO REACH COMMON RATIO ■ 21 
NljiRFR OF RESPOnCF FaCTDR^ PFR _5ET '• 22 _ 

CO wj^mT RATIO _• *_7 3 NB32 Q9 69 
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table 8 (Cont.) 


NEW RES FAC 


^r^ONSE FACTORS' 


HOUR 

X 

Y 

~ I 

0 

.7251569*15 

.0003756761 

.8374457390 

1 

—.6562617794 

.0036660957 

-.5326726911 

2^ 

-.0101950074" 

.0055020236"" ’ 

-.1759137619 

3 

-.0074798671 

.0050052263 

-.0596083779 “ 

4 

-.0055072146 

.0039969156 

-.0211851230 

5 

-.0040565345 

".0030453759* 

-.0082261792 

6 

-.0029085459 

.0022765110 

0036634233 

7 

-.0022019155 ■” 

.0016879889 

—• 0C1 9206029 

e 

-.C016??3O73 

.0012472019 

-.0011622680 

9 

-.0011954207 

.0009200953 

.-.0007742684 

lo ' 

-.0008306205 

.0006783189 

-.0005438347 

n 

-.0006*90160 

.0004999255 

-.0003920536 * 

12 

—.00047*21 60 

• COC3603998 

— • 0C02860651 

13 

‘-.00035*3651 

.0002714615 

-.0C02098687 

14" 

-.0002596342 

.0002000259 

-• 0C01 543415 

15 

-.0001913070 

.0001473870 

-.0001136275 

16 

-.0001409613 

•0001086000 

-.0000836932 

' I?" 

-.0001038649 

•OO0C8OO2O2 

-.0000616578 

IP 

-.0000765311" 

.0000589616 

-.0000454282 

l<j ‘ 

-.0000563907 

" .0000434449 

-.0000334719 

2C 

-.000041 *505 

.0000320116 

-.0000244629 

21 

-.0000306158 

• 0000235672 

-.0000181723 


" NUM»f R OF HOl'kS RFOUIRFD TO REACH C ON HON RAT 10 • 21 

'^•M«rir-cr-7E5POK , CE ‘FACToirr'P r FirsET ■ 22 

'to'nHON patio' ■ ' ‘.7368330005 "" 
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TABLE 8 (Cont.) 


ROOTS 


_ ROOT « - OLD 


NEW 


1 0.30539401103140 

2 1.12468525501743 

3 9.41824542784383 

-A - — 12.27360260558755 

5 — -33.85678581614047 


0.30539400613318 

1,12468526089836 

9,41824542612801 

12.27380260257945 

33.85678582246646 


b (ROOTS) 


FOR ROOTS ft OLD 


NEW 


1 -3.4314850E-7 

2 -2.9825105E-7 

3 -5.2242660E-8 

A 7.2314770E-& 

— 5 - 6.6737980E-8 


2.93000E-12 
-1.36000E-12 
-7.50152E-9 
-1 .06420E-10 
3.69000E-12 
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TABLE 9 


DESCRIPTION Of' CONSTRUCTION 


LAYER 

_ THICKNESS. 

conductivity 

DENSITY 

. . SPECIFIC. HEAT 

.. RESISTANCE _ 


NUMBER 

FT 

BTU PER t HR I C FT1 <F > 

LB PER CU 

FT BTU PER ( LB M F ) 

f HR) (SO FTHF1 PER 

BTU 

1 . _ 

*0417 

.830 

55.0 

.4000 

0,0000 

STONE 

2 

• 0313 

.110 

70.0 

• 4000 

0.0000 

FELT 

3 

.1670 

.025 

2.0 

.2000 

0.0000 

_ CELLULAR .CLASS 

4 

_ .3330 

.100 

40.0 

.2000 

0.0000 

_ CONCRETE 

5 

.0030 

26.000 

480.0 

• 1000 

0.0000 

metal pan 

6 

o.oooo* 

0.000 

0.0 

0.0000 

.0100" 

AIR SPACE 

7 

" .0625' 

• 420 

“100.0* 

• 2000 

0.0000 

GYPSUM BOARD 

8 

0.0000' 

' 0.000 

0.0 " 

0.0000 

.9100 

AIR SPACE 

9 

.0625" 

.035 

“' 30.0“ 

• 2000 

0.0000 

" ACOUSTICAL TILE 

10 

0.0000' 

0.000 

0.0 

0.0000 

• 6850 

INSIDE AIR 


OLD RESFAC 


RES PON SE ~ f AC TORS"' 


HOI'Q 

X 

Y 

2 

0 

.757R139001 

.0027676836 

.845 5636399 

1 

-.3264110204 

.0257121522 

-.3394234703 

? 

-.0551369293 

.0378060569 

-.1224280305 * 

3 

-.0A06A6719A 

.0340939470 

-.0561552755 

A 

-.0303621399 

.0272297798 

-.032 1965209 

5 

-.0227925210 

.0209383902 

-.021 4183746 

6 

-.0171A1RO60 

.0150870616 

-.0153359574 

7 

-.0129010796 

.0119950160 

-.01 1 3254 352 “"* 

8 

-.0097119210 

.0090A13992 

-.0004651492 

q 

— • 00 7 3 1 1 0 2 7 A 

♦0040 100783 

-.0063561928 

10 

-. 005305061 A 

.0051281546 

-.004 7308134 

ll . 

-.0GA1AAR055 

.0030612698 

-.0035981 829 

12 

-.0031206732“ 

.0029072631 

"-.002 7067421” 

13 

-.002349 5961 

'.0021889357' 

- • 002 C3934 3 5 

1 A 

-.0017690A32 

.0016400851 

-.001 5354205 

15 

-.0013319373 

".'0012408677“ 

-.001 1560313 






KIMPRFP pp Mnu»S RFOtJIfifn Tn PFACH common RATIO • 15 



N 1 1 N • F 0 pc PFSPONCF FATTOOS PFP 

SFT • 16 



cn**“nn RATIO • * 752 °1 4Q4 5 2 
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TABLE 9 (Cont.) 


NEW RESFAC 


^50QN5E>AetO*if 


HOUR 

o 

j ■ - 

X 

•7570138940 


2"“ "" 

—.3264 1 10153 
-.0551369705 


3 

4 

-.04084671 89 
-.0303671396 


5 

-.0227925708 


6 

-.0171414959 


7 

-.0179010796 

— 

n 

—.009711 921 9 


9 

-.0073118774 


10 

-.0055050615 

■ — — • 

n 

-.0041448056 

— 


H-.0031 206733 

— 

13" 

-.0023495962 

. . ... < 

14 

—•GO 17 6904 33 

1 

15 

-.0013319373 

__ a 


• 002 76770 u] 
•0297121323 
.0378860570“ 
^ 0340939471 ' 


_ JJ0435636374 
^-.3394234430" 
^.1724280295 
0561552749“ 


.0272297800 -•0371965207 

.0209303903 -.0214181745 

.0158870614 -.C 153359574" 

• 0119950162 -.0113254352 " 


• 00904 13994 " 


.0051701348" 
• 00386 12700 
.0029072632 
.0021589358' 
.0016400851 " 
,00124086771 


-.0084651493" 

“>0063561929 

-.0047804135 

-.0035981530 


-.0020393435 

^.0015354206 

_^,00115603J4^ 


— NUHBFR Of HOURS FEOUXPED "Vo'ltFACH Com HON •ATItP* 15 * 

NUIBfP PESPbNCE f ACTORS ' pTp'SFT"* 16 

COW*PN PATIO »" "*V 7 5291404 6 3~ 


ROOTS 


ROOT-# OLD - - — NEK* - - 


1 0.-28380420751416 0 .28380420343971 - 

2 1.27316441620625 1.27316442200970 - 

3 6 .-089163X2 121 456 6.089163722X0650 

4 16.5953121X201106 16.59531217589722-- 

— 5 31^14535162877291 31,14535163507094-- 


B(rOOTS) — • 


^OR-ROOT-S-» OLD - _ NEW 


1 -4.681876E-8 1.800E-12 

2 »4.410708E-8 -3.000E-14 - 

3 : 1.978496E-8 --1.238E-11 

4 -*5.435033E«8 -8.071E-1-1 

5 3.95i397E=8 _ . . . 1.320E-12 
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TABLE 10 


1 ..1*1 DE S CRI P T I ON "OF . CONST RUCTION. 


L AYER TH1 CKNF SS ■ C ONDUCTIVITY DENSITY SEEdF-IC-HEAT RESISTANCE 

-X*U W 8CB FT 0TU PER (HR) (FT) IF) LB PER CU.FT BTU PER (LB) (F) (HR) (SU FT) (F) PER- BTU 

1 *0050 26.000 *80 • 0 *1000 - 0.0000 METAL SIDING 

Zi .0*20 .025 2.0 .2000 0.0000 INSULATION-PANELS — 

3 .0«33 .*00 116.0 .2000 0.0000 STUCCO 

4 .0313 .110 .70.0 .4000 0.0000 FELT 

5 .0420 .067 34.0 .2900 0.0000 PLYWOOD .SHEATHING 

6 0.0000 0.000 0.0 0.0000 .9100 AIR SPACE 

7 .0625 .420 100.0 .2000 0.0000 GYPSUM BOARD 

$ .0104 .070 37.0 .6000 0.0000 PANELING 

9 0.0000 0.000 0.0 0.0000 .6050., INSIDE AIR 


OLD RE S FAC 


HftESPONSJ? /ACTORS 


HOUR 

X 

Y 

Z 

0 

1.8063060277 

• 0000004^10 

.4 01 1426131 

1 

-1.7737510020 

.0001551499 

-.151 2127477 

? 

-.0003253801 

.0012842006 

-.0294006032 

3 

-.00553fr«512 

.0028710285 

— • 020 194204 6 

u 

-.0043045688 

.0036685829 

-.0152586280 

5 

-.0035601171 

.0042707124 

-.012 2500417 

6 

-.0030526487 

•0043131968 

-.0102674603 

7 

-.0026786892 

.0041663052 

-.008 66938 7 1 

p 

-.0023856600 

.0039289601 

-.0078185899 

Q 

-.002X447051 

.0036553284 

-.0069833746 

10 

-.001O3O347R 

.0033741021 

-.0062892735 

11 

-.0017^99490 

•0031003295 

-.005 6933503 

12 

-.0016006431 

• 002 64 0 88 3 4 

'-.0051701766 

13 

-.0014576938 

.0025988771 

-.0047041173 

1 * 

-.0013285809 

.0023751*41 

-.0042850716 

15 

-.0012114942 

.0021693653 

- • 003 9061 i 6 F 

16 

-.0011050518 

.0019807441 

-.003 5621987 

17 

-.0010081407 

• 001 60P 11 72 

-.0032494010 

lfi 

-.0009108274 

.0016503210 

-.0029645329^" 

10 

-.0008303047 

• 001506 1780 

— .002 70 489 3 5 ~ 

20 

-.0007658611 

.0013745590 

"-.0024681341 

21 

-.0006988606 

.0012544076 

-.0022521756 

22 

—.0006377307 

,0011447383 

‘-. 002 0551558 

23 

-.0005810520 

.0010446462 

"-.0018753946 

24 

-.0005310563 

.0009532099 

-.001 711369B 


RUK ftFP p c HOURS RFQHIR£D_TO PEACH COMMON "RATIO ■ 

nTjm'bj Q PS PONCE' e ACTOR 5 PF R_S E T • 25 

COM MON PATIO V ’ .9 12 54 AO 730* * 
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Pfc<P~OK *r FACTO* $ 



NUK0ER OF HGUPS RECUIRED TG *6 KiTCO^PN '» ATIO • 24 

~7CMRER“DF”PESi>CWE"PACTOPS PTR SET ■' 25 

*~C 0 KHCN "PiTICr* "”.'91254 6 9711 
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TABLE 10 (Cont.) 


ROOTS 


- ROOT—# — : OLD-- 



-0.-09151352592139 - 

—0.59670697196906 - 

-1.87826524925806 
—4.46784784597435 
-5.21288150152188 
-6,84838752278688 
—9.19016020960470 
15.55378495246231 
21,54078931780691 
25,48150421908758 
26,15232440552245 
36,39218806754798 


- b (ROOTS) - 


FOR- ROOTS-# 


-OLD 


1 2.41 13701E-7 

-• — 2 2.6229038E-7 

3 1*005886 1E“7-- 

4 -1.3448400E-9 

5 8.6910700E-9 

_ 6 1 .8604440E-8 

-7 -1.0867190E-8 

. .. -a -4.9799730E-8 

~ — 2.9246420E-8 

- -10 5.6764700E-9 

11 - 7.1909600E-9 

12 - - — 3.9220780E-8 


NEW 

- 0.09151352798700 
0.59670696614259 
1.87826525339801 

- 4.46784784772299 
5.21288150553909 

6,84838751809579 

- 9.19016020792799 
-- 15.55378495460405 
- 21.54078932206482 

25.48150421445916 
-- 26.15232441059038 
- -36.39218806540566 


NEW 

1.50000E-13 
-2.40000E-13 
-9,3l998E-9... 
4.2862QE-9 
1.92000E-12 
-2.90600E-11 
7.27528E-9 
1.23000E-12 
-2.2931 3E-9 
9 .40000E-13 
2.68070E-10 
-5.25500E-11 
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APPENDIX C 


BISECTION METHOD FOR LOCATING ROOTS 

Given a continuous function f on the interval [a, b] , where f ( a ) 
and f(b) have opposite signs. To find a solution to f ( x ) = 0.0: 

1) set a j = a and b j = b 

2) set i = 1 

3) set = Ma.+b.) 

4) if ABS( f ( p.. )<1 . OE-8 go to 6, else go to 5 

5) if f(p i -)*f(a i )>0,set a. + 1 = p.. and b i+i =b. ; 
if f(P i )*f(a i )<0,set a i + 1 =a. and b 1 + 1 = p i ; 
add 1 to i and go to 3 

6) end 
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APPENDIX D 


SECANT METHOD FOR LOCATING ROOTS 


Given a continous function f on the interval [ 
and Pi are initial root approximations such that Po 
a solution to f(x) = 0.0: 

1) set i = 2 


2) set p n> 


f (Pi.i ), (p i . 1 -P i - 2 ) 

f( Pi-l)- f (Pl-2> 


Po 9 

\ 


P 

Pi 


il 


, where p 0 
To find 


3) if ABS(f(p^)) <_1 . 0 E - 8 go to 4, else add 1 to i and go to 2 

4) end 
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